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ABSTRACT 

We report Mopra (ATNF), Anglo- Australian Telescope, and Atacama Submillimeter Tele- 
scope Experiment observations of a molecular clump in Carina, BYF73 = G286.21+0.17, which 
give evidence of large-scale gravitational infall in the dense gas. From the millimetre and far- 
infrared data, the clump has mass ~ 2xlO^M0, luminosity ~ 2-3xlO^L0, and diameter ~0.9 
pc. From radiative transfer modelling, we derive a mass infall rate ~ 3.4x10"^ M0yr~^. If 
confirmed, this rate for gravitational infall in a molecular core or clump may be the highest yet 
seen. The near-infrared i^^-band imaging shows an adjacent compact HII region and IR cluster 
surrounded by a shell-like photodissociation region showing H2 emission. At the molecular infall 
peak, the K imaging also reveals a deeply embedded group of stars with associated H2 emission. 
The combination of these features is very unusual and we suggest they indicate the ongoing for- 
mation of a massive star cluster. We discuss the implications of these data for competing theories 
of massive star formation. 



Subject headings: astrochemistry — infrared: ISM 
— radio lines: ISM — stars: formation 



ISM: kinematics and dynamics — ISM: molecules 



1. Introduction 



2002 



Many details of massive star format ion in dense 

molecular clouds are still unclear ( Churchwell 

2002 ) , despite much recen t progress (e.g., Sridhara n et al.l 
iFuller et all l2005l : iLongmore etal. 2007). 



For example, it is still debated whether massive 
stars can form by a scaled-up version of the ac- 
creti on thought to occur w ith low-mass protostars 
(e.g. iMcKee fc fanl 120031 ). or rather form by col- 



1 School of Physics A28, University of Sydney, NSW 
2006, Austraha 

^Astronomy Department, University of Florida, 
Gainesville, FL 32611, USA 

^Department of Physical Science, Osaka Prefecture Uni- 
versity, 1-1 Gakuen-cho, Sakai, Osaka 599-8531, Japan 

^Center for Astronomy, Ibaraki University, 2-1-1 
Bunkyo, Mito, Ibaraki 310-8512, Japan 

^Anglo-Australian Observatory, PO Box 296, Epping, 
NSW 1710, Austraha 

^Department of Astrophysics, Nagoya University, Furo- 
cho, Chikusa-ku, Nagoya 464-8602, Japan 



lective processes in a clustered environment (e.g. 



Bonnell et al. 20031 ). Consequently, examples of 



massive star formation showing evidence of either 
behaviour can be informative to this debate, espe- 
cially since there are still relatively few examples 
known of true massive protostars. 

As part of the Census of High- and Medium- 
mass Protostars (CHaMP, §2.ip . we identified the 
massive dense clumpE] G286.21+0.17 as showing 
striking evidence of large-scale gravitational infall, 
which we report here. This source (hereafter re- 
ferred to as BYF73, from the master CHaMP 

I2OO9 ) has been in- 



sour ce list; Barnes et al 



elude d in some prev i ous s ur yevs (iBronfman et al" 



1996: Dutra et al 



2OO3I : iFaund ez et all l2004l : 



Yonekura et al.l 12005 ) and is an Infrared Astro- 
nomical Satellite (IRAS) point source, but has 



^Here we use the IWilh ams et"aD (120001 ) terms for " core 
(that part of a molecular cloud which will collapse to form 
an individual star or binary) and "clump" (which will col- 
lapse and fragment, via many cores, to form a star cluster). 
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not previously been shown to be remarkable. The 
precise location is {l,b) = (286°208, +0°169) or 
{a,6)j2ooo = (10^38^3212, -58°19'12''), about 
1°5 northwest of rj Carinae and 12' north of 
the rim of the 15'-diameter HII region/bubble 
NGC3324/IC2599, at an assumed distance of 2.5 
kpc. 

2. Observations 
2.1. Survey Strategy 

The motivation for CHaMP is to make a com- 
plete and unbiased census of higher-mass star for- 
mation at many different wavelengths over a large 
portion of the Milky Way (Barnes et al. 2006) . 
in order to systematically characterise the pro- 
cesses in massive star formation in a uniform 
way. The first step was to identify 209 dense 
clumps from C^^O and H CO+ maps made with 
the 4m Nanten telescope (Yonekura et aL 20051 . 
[2009) of a 20° X 6° region of the Galactic Plane in 
Vela, Carina, and Centaurus (specifically 300° > 
/ > 280° and -4° < 6 < +2°). A higher- 
resolution follow-up campaign was then begun to 
map these clumps in a number of 3-millimetre 
wavelength (3mm) molecular transitions with the 
22m-diameter Mopra dis h of the Australia Tele- 
scope National Facilit}|l (|BarneseEaD[20^ The 
Mopra antenna's performance has been described 
by Ladd et al. (2005). Since that study, an on- 
the-fiy (OTF) mapping capability has been imple- 
mented in the control software (in 2004), new 3mm 
MMIC receivers were installed (in 2005) which 
were at least as sensitive as the previous SIS mix- 
ers and much more efficient to operate, and the 
MOPS wideband digital filterbank was commis- 



sioned (in 2006; IWilson et aP l2006h . These in- 



novations, when combined with the Nanten maps 
as finder charts, makes an ambitious survey like 
CHaMP possible. 

Mopra's OTF mapping mode has been de- 
scribed by T. Wong (2005, unpublished), which 
we briefly summarise here. The telescope is driven 
in a raster pattern (which can be in any of the 



^The Mopra telescope is part of the Australia Telescope 
which is funded by the Commonwealth of Australia for op- 
eration as a National Facility managed by CSIRO. The Uni- 
versity of New South Wales Digital Filter Bank used for the 
observations with the Mopra telescope was provided with 
support from the Australian Research Council. 



6, a, or 5 directions) at a rate such that the data 
dump interval (usually every 2 s) from the correla- 
tor to mass storage is consistent with Nyquist- (or 
better) sampling of the sky, given the telescope 
beam and observing frequency. At 90 GHz this 
drive rate across the sky equates to approximately 
6'' sec~^. Each raster row is then offset by a simi- 
lar amount (i.e., 12'' at 90 GHz) from the previous 
row, until a square map with a size of the user's 
choosing is built up. The user also selects whether 
a reference position (which can be specified in ei- 
ther relative or absolute coordinates) is observed 
at the beginning of each row, or only once every 
2 rows. Additionally, the user can choose from 
which corner of the square map the raster pattern 
is begun, i.e. the NE, NW, SE, or SW (in the re- 
spective coordinate system being used). Finally 
the frequency of hot-cold load measurements of 
Tsys needs to be specified; this is typically ev- 
ery 10-30 min, depending on the stability of sky 
conditions. In the 2007 season, however, a noise- 
diode calibration system was introduced into the 
data stream, effectively giving continuous Tsys 
measurements and making separate hot-cold load 
scans somewhat redundant. With the addition 
of a calibration spectrum of a known source such 
as Orion-KL, skydip measurements of the atmo- 
spheric opacity were not found to be necessary. 

In this way a typical 5'x5' map can be built up 
over a period of about 70 min at 90 GHz. In or- 
der to minimise rastering artifacts, however, a sec- 
ond map is usually made of the same field, but in 
an orthogonal rastering direction. Including time 
(~ 10 min) for pointing checks between each map, 
such a b'xb' field is "complete" in about 2.5 hr. 
Further rasters can be done of the same field, and 
this not only improves the S/N in the usual way, 
but under variable sky conditions will also min- 
imise noise variations across a map, which might 
otherwise give erratic sensitivity coverage of the 
user's field. After just 2 raster maps, however, the 
noise variations are usually acceptable (^20%) in 
all but the worst conditions. 

The MOPS backend can be employed in either 
"broadband" or "zoom" mode. With the former, 
the full 8 GHz available bandwidth is observed 
with 65536 125-kHz-wide channels in each polar- 
isation, corresponding to a velocity resolution of 
0.45 km s~^ at 90 GHz. In contrast the latter al- 
lows up to sixteen independently selectable 138- 
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MHz-wide "zoom IFs" to be observed simultane- 
ously from within the filterbank's 8 GHz total in- 
stantaneous bandwidth. Each zoom mode is cor- 
related with 4096 channels in each of two orthogo- 
nal polarisations, resulting in a spectral resolution 
of 33 kHz, or 0.11 km s"^ at 90 GHz. In the 2005- 
07 austral winter seasons we mapped the brightest 
118 Nanten clumps, simultaneously covering many 
spectral lines in the 85-93 GHz range, among them 
the J=1^0 transitions of HCO+, HCN, N2H+, 
H^^CO+, and H^^CN. At these frequencies, Mopra 
has a beam FWHM of 36'', an inner error beam 
which extends to ~80'', and a coupling efficien cy 
of 0.64 to sources of this size (|Ladd et al.ll2005h . 

While CHaMP's 3mm molecular maps reveal 
the location of dense gas, complementary near- 
IR imaging of the same clumps can show where 
star formation has evolved further. By compil- 
ing these statistics uniformly we will be in an 
excellent position to identify demographic trends 
in the massive star formation process. Thus, an 
equally important part of CHaMP is a near-IR sur- 
vey of the Nanten c lumps using the IRIS 2 imager 
( Tinnev et al ] l2004h on the Anglo- Australian Tele- 
scope (AAT). With this instrument we have begun 
acquiring images of each clump in i^-continuum, 
Brackett-7 (a recombination line tracing HII re- 
gions), and H2 ^=1^0 S{1) & ^=2^1 S{1) (vi- 
brational quadrupole lines tracing molecular gas 
heated to a few 1000 K) to delineate the relation- 
ship between formed and forming massive stars, 
and report here results of such imaging. 

A third major component of CHaMP will be 
a deep imaging survey of 1.2mm continuum and 
spectral-line emission with Atacama Submillime- 
ter Telescope Experiment (ASTE), the 10m sub- 
millimeter telescope of the Nobeya ma Radio Ob- 
serva torJl at Pa mpa la B ola in Chile feohno et al 



2004t [Ezawa et al]l2004f ). The 1mm continuum 



is important in characterising the spectral energy 
distributions (SEDs) of embedded protostars as 
well as starless cores or clumps, and in correlating 
this with the phenomenology seen in spectral lines 



^The ASTE project is led by Nobeyama Radio Observa- 
tory (NRO), a branch of National Astronomical Observa- 
tory of Japan (NAOJ), in collaboration with the University 
of Chile and Japanese institutes including the University 
of Tokyo, Nagoya University, Osaka Prefecture University, 
Ibaraki University, Kobe University, and Hokkaido Univer- 
sity. 



and at other wavelengths. We report here as well 
some of the first data from this facility, namely 
HCO+ and R^^CO^ J=4^3 spectra, confirming 
the evidence of infall from our Mopra data. 

2.2. Observational Details and Data Re- 
duction 

The evidence for infall in BYF73 was first seen 
in the Mopra HCO+ and }1^^CO^J=1^0 data, 
presented in Figures [THU These were obtained on 
2006 Oct 27-29 and 2007 Sep 5-9, when observ- 
ing conditions were good {Tgys ~ 300 K or better). 
The images were formed by coadding 5' x 5' OTP 
maps which abut each other to cover larger areas. 
The reference position used for sky-subtraction 
during all mapping was = (285.7,-0.3), which 
shows no emission in the Nanten CO map. Each 
5' X 5' area was scanned twice or three times in 
each of / and b in order to minimise rastering ar- 
tifacts and noise variations. The raw OTF data 
were processed with th e Livedata-Gridzilla pack- 
age ( Barnes et al.ll200lh by bandpass division and 
baseline subtraction. The 2s-long OTF samples 
were then regridded onto a regular grid of 12^' 
pixels, where the samples were weighted by T^l, 
before averaging them into each gridded pixel. 
Weighting by the rms~^ of the spectra was not an 
option provided by Gridzilla; however as described 
above, since 2007 the continuously-measured Tsys 
has effectively given the same information for each 
2 s sample. For all Mopra maps in Figures [T]-[7j the 
effective telescope HPBW has been smoothed at 
the gridding stage to 40'' from the intrinsic 36'', 
in order to reduce noise artifacts. The resulting 
spectral line data cubes have low but, due to vari- 
ations in weather and coverage, somewhat variable 
rms noise levels, ranging from 0.17 K in the south- 
ern portions of Figure [2] up to 0.22 K in the north, 
per 0.11 kms~^ channel; the average across the 
map is 0.2 K per channel. Although the pointing 
(checked on the SiO maser source R Carinae every 
hour or two) was typically good to 10" or better 
(<1 pixel on the scale of our maps), because of the 
simultaneity of the spectral line mapping afforded 
by MOPS, the relative registration of features be- 
tween these lines is perfect. 

Observations of HCO+ and H^3CO+J=4^3 
were made using ASTE on 2006 Dec 1-2, when 
the typical system temperature (double-sideband) 
ranged from 220 K to 580 K at 345 GHz, includ- 
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286.25 286.2 286.15 286.25 286.2 286.15 286.25 286.2 286.15 286.25 286.2 286.15 286.25 286.2 286.15 

GLON (degrees) GLON (degrees) GLON (degrees) GLON (degrees) GLON (degrees) 

Fig. 1. — Mopra channel maps of HCO+ J=1^0 emission from an inner 7'-field of BYF73 = G286.21+0.17, 
shown on the T| scale; = T|/0.64. Each panel is actually the mean of four 0.113 kms~^ wide channels, 
hence the panel spacing of 0.453 kms~^. The grey scale is linear from OK to the peak at 2.19 K, while the 
contour levels are from 3 to 9 times the rms noise level for each panel, which is 0.092 K. At a distance of 
2.5 kpc, the scale is 40'' = 0.485 pc or 0°.02 = 0.873 pc or Ipc = 0°.0229 = 82'!5. 



ing the atmosphere. The half-power beamwidth 
of the telescope is 22'' at 345 GHz and the front 
end is a 4 K cooled SIS mixer receiv er; at this fre- 
quency the beam efficiency is 0.65 (iKohndliooHl : 



Kohno et al. 20081 ). We used a digital correlator 



with a bandwidth of 128 MHz and 1024 channels 
( Sorai et "al]|2000l ). The effective spectral resolu- 
tion was 151.25 kHz, corresponding to a velocity 
resolution of 0.13 km s~^ at 345 GHz. The data 
were obtained in position switched mode centred 
on {l,b) = (286°2071, +0°1692); the "off" position 
used for sky-subtraction was (286°.0407, +0.4612), 
which is also devoid of Nanten CO emission. Ob- 
servations were made remotely from an ASTE op- 
eration room in San Pedro de Atacama, Chile, us- 



ing a network observ ation system, N- COSM O S3, 
developed by NAOJ ([Kamazaki et all [20051 ). For 
HCO+ and H^^CO^, the total integration times of 
the spectra were 180 s and 620 s, and the rms noise 
levels 0.57 and 0.21 K per channel, respectively. 
The intensity was calibrated by using a room- 
temperature chopper wheel. The absolute inten- 
sity was calibrated by observing Orion-KL and as- 
suming T*(HCO+ ) = 4 7 K and T*(Hi3C0+) = 
2K ([Schilke et anil997h . The pointing accuracy 
was measured to be reliable within 5" as checked 
by optical observations of a star with a CCD cam- 
era attached to the telescope, as well as by CO 
J=3^2 observations of IRC+10216. 

IRIS 2 observations were made at the A AT in 
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service mode on 2006 May 13 in the iC-band 
continuum (2.25-2.29 /im) and 3 spectral line 
filters as described above, Br-7 and H2 
'^=1^0 & v=2^1. For each filter, nine 60s im- 
ages (dithered by 1') were obtained of the instru- 
ment's I'.lxl'.l field under 0'!9 seeing, and reduced 
using the ORAC-DR data reduction pipeline 
( Cavana gh et al.^'2003^. However the conditions 
during the observations were non-photometric, 
there being a fair amount of bushfire haze present. 
Subsequent image processing was performed with 
the IRAfH package. The images for each field were 
registered using astrometry derived from Super- 
I-band images; we estimate the result- 
ing rms positional accuracy in the IRIS2 images 
to be < 0'!3. Next, we linearly scaled the spectral- 
line images to the same relative brightness scale as 
the K-band continuum by matching the integrated 
fluxes of several stars in each filter, assuming they 
were of similar colour. We then subtracted the 
continuum from the spectral-line images before 
transforming each image to Galactic coordinates. 
Finally, a three-colour image (shown in Fig.[7K) 
was formed from the continuum-subtracted Br-7 
and H2 images: Br-7 is shown as red, and H2 5'(1) 
'^=1^0 & v=2^1 are shown as green & blue, 
respectively. 

Long-slit spectroscopy with IRIS2 was obtained 
on 2007 Oct 18. The 7!7 long slit was set to a po- 
sition angle of 131°2, with the stellar cluster and 
nebulosity of BYF73 spanning most of one half of 
the slit. Four exposures of 300s were obtained in 
the i^-band, with BYF73 nodded by 3!8 along the 
slit between each exposure. Similar nodded expo- 
sures of the nearby AOV star HD 95534 were ob- 
tained to assist with telluric correction. All frames 
were flatfielded using quartz lamp exposures, then 
nodded pairs were subtracted to remove sky emis- 
sion. After two-dimensional wavelength calibra- 
tion and straightening with Xe lamp exposures, 
the "off" beam data were inverted, aligned, and 



HCO+ and H^^CG^ integrated intensities 



+0.20' 




+0.05' 



286.30° 286.25° 286.20° 286.15° 

Galactic Longitude 

Fig. 2. — {Image) Wider ll'-field of Mopra 
HCO+J=1^0 integrated intensity from BYF73, 
on the scale as given by the colour bar. 
The integration is over the range -23.20 to - 
16.63 kms~^ or 58 channels, yielding an rms noise 
level 0.16 K kms~^: hence the widespread low- 
level emission above ~0.5K kms~^ is real. {Con- 
tours) Mopra H-^^CO+J=1^0 integrated inten- 
sity in levels are {grey) -0.5, -0.35, {black) 
0.35, 0.5, 0.7, 0.9, and 1.1 K kms-^ The integra- 
tion is from —21.94 to —17.86 kms~^, giving an 
rms noise level 0.12 K kms~^. The smoothed Mo- 
pra HPBW for both datasets (40'') is shown for 
reference in the lower-left corner. 



co-added to the "on" beam data. Each spectral 
row of the data was divided by an extracted spec- 
trum of HD 95534 (from which intrinsic Br-7 ab- 
sorption had been removed), then multiplied by a 
blackbody spectrum of Tg// = 9520 K. 



^IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under coop- 
erative agreement with the U.S. National Science Founda- 
tion. 

^This research has made use of data obtained from the Su- 
perCOSMOS Science Archive, prepared and hosted by the 
Wide Field Astronomy Unit, Institute for Astronomy, Uni- 
versity of Edinburgh, which is funded by the UK Science 
and Technology Facilities Council. 



3. Analysis and Discussion 
3.1. Mopra Maps 

The Mopra HCO+ J=1^0 maps, being of high 
signal-to-noise ratio, reveal a number of interest- 
ing features which we describe here. Figure [T] 
shows the HCO+ emission from BYF73 across 
its full velocity range, where we have averaged 



5 



(a) (b) 




Galactic Longitude Galactic Longitude 

Fig. 3. — Moment images of BYF73 from the Mopra HCO+ data with telescope beam as shown in the corners, 
and overlaid with HCO+ integrated intensity contours at 1.35, 1.9, 2.4, 3.0, 4.0, 5.0, and 6.0 K kms~^ (a = 
0.16K kms-^). Ah moments for HCO+ were calculated over the same velocity range as in Fig.^image). 
(a) First moment (intensity- weighted mean velocity field), (b) Second moment (velocity dispersion). 



four velocity channels into each displayed panel 
for ease of viewing. However all analysis below 
rests upon the full-resolution data. At the central 
velocities (-22 to -18kms~^) the emission is quite 
widespread; redward of the line centre (-20 to - 
18kms~^) this extended emission is quite clumpy, 
while to the blue (-22 to -20kms~^) the emission 
is strongly centrally concentrated. At both the 
reddest and bluest velocities, the emission is fairly 
centrally concentrated; in particular there is lit- 
tle obvious evidence for an extended, high- velocity 
outflow which would tend to have emission well 
away from the centre at the highest relative ve- 
locities. In fact to a casual inspection there seems 
to be little systematic kinematic structure to these 
channel maps at all, and it is only upon inspection 
of the spectra that the infall profiles are revealed. 

This is also reflected in the integrated intensity 
image of Figure O where the extended envelope of 
BYF73 shows very little evidence of being struc- 
turally disturbed by (for example) its proximity 
on the sky to r] Carinae or to NGC 3324. The 
only morphological feature of note in the envelope 
is a bay to the NW, which as we shall see is in- 
trinsic to the source. The inner 2' of BYF73 also 
appears fairly bland: this area is elongated some- 
what in the NW-SE direction in both the HCO+ 
and H^^CO^ emission, and there is a small but 



significant offset in the peak positions of these 
two molecules, with the H^^CO^ emission centred 
slightly to the northwest of the brightest HCO+. 

In Figure [3] we give the higher-order HCO+ 
moment images, overlaid by the moment-0 con- 
tours from Figure [21 In contrast to the latter, 
the intensity-weighted mean velocity (Fig.[3]2) re- 
veals a striking velocity gradient across the clump; 
the axis of this gradient is rotated ~ 30° anticlock- 
wise from the long axis of the clump, as seen in 
the moment-0 contours. The spectral line is most 
strongly blue-shifted to the north and east of the 
peak HCO+ emission, reaching its minimum value 
^0^5 east of the peak. This blueshift gradually 
changes to a redshift to the western side of the 
clump, reaching its maximum value ~2' west of 
the peak, inside the bay of the envelope where the 
emission is weaker. The weaker emission which 
wraps around the western side of the bay (the 
"western arm") is also reshifted with respect to 
the clump. As a whole the clump's velocity is sig- 
nificantly blueshifted with respect to the H^^CO^ 
line centre (see below), which is approximately at 
the green colour in this image. 

The IICO+ velocity dispersion a (where the 
line FWHM = 2.355cr) is shown in Figure 
[Note that because the line shape is strongly non- 
Gaussian, one should not confuse this moment-2 
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measurement with an actual linewidth; neverthe- 
less its variation does indicate true changes in the 
line profile across the source.] Here again we see a 
strong gradient in this parameter: the bulk of the 
clump, and to its north and east down to an inte- 
grated intensity level ~2.0Kkms~^, has a large 
cr ~ 1.5-1. 8kms~^. Below this intensity to the 
west and SW, a drops to 1.1-1.4 km s~^, reaching 
minimum values <lkms~^ exactly where the ve- 
locity field is most redshifted. One therefore sus- 
pects that these redshifted features are due mostly 
to individual, narrower-line substructures in the 
envelope, and that apart from these features, the 
clump's overall blueshift with respect to the opti- 
cally thin line centre is even more complete. 

These images can be compared to the equiv- 
alent moment images of the H^^CO^ cube (not 
shown here, but see Fig.Hf) for a spectrum). From 
Gaussian fits to this cube, the peak T\ ~0.3- 
0.4 K, and the line is centred near Vlsr ^ ~20 to - 
19.5 km s~^ across most of the emission, but shifts 
to -19 to -18.5 km s~^ along the clump's SW edge. 
The line FWHM varies from ~ l-3kms~^ toward 
the SE end of the clump, rising to ~2-5kms~^ 
toward its NW end. Such fits to the H^^CO+ 
have rms residuals ~0.17K per 0.12 km s~^ chan- 
nel, with typical uncertainties ~ 0.3 and 0.7kms~^ 
to the Vlsr and linewidth, respectively. 

3.2. Distance Determinations 



Typical Galactic r otation curves ( Burton fc Gordonl 
19781 : Icieme"^ Il985h and standard values of 
and Go (8.5 kpc/220 kms-\ lAU 1978 values; 
8kpc/200 kms~ \ Merrifield 1992; 8.4kpc/254 kms-\ 
[lleid et al.l l2009l ) indicate that the central Vlsr 
= -19.7 kms~^ for the molecular clump (see §3.4p 
is formally forbidden at the longitude of BYF73 
(meaning that this velocity is inconsistent with 
such rotation curves for objects at any distance 
along this line of sight). In fact the minimum 
allowed velocity for the tangent point at this 
longitude is ~ — lOkms"^, which is ^2a more 
positive than Vbyfys (where a = 4 kms~-^ is 
the cloud-to-cloud velo city dispersion of GMCs; 
Burton fc Gordonll 19781 ) . In spite of this disparity, 



any other location for BYF73 is even less kinemat- 
ically favoured than the tangent-point distance. 
For example, BYF73 and the r] Car GMC may 
be sharing in a non-circular streaming motion of 
the order of 10 kms~^ associated with this part 



of the Carina Arm. In any case, at such tangent 
points small uncertainties in the rotation curve or 
the values of Rq and 6o can translate into large 
line-of-sight distance uncertainties, up to 50% or 
more. Therefore while a tange nt-point distance 
of i?ocos(/)=2.35±1.5kpc fusing iReid et al.ll2009h 
is favoured with the kinematic method, a more 
robust determination is preferred for the analysis 
in the sections following, especially in light of the 
large power with which the distance to BYF73 
enters some of the formulae below. 

Fortunately, a number of studies have yielded 
distances to the massive clusters in and near 
the 77 Car GMC ( e.g., see the sum mary by 
Yonekura et"alll2005h . and to NGC 3324 ([Havnes et al 



1978[ ). These range from 2.2-2.8 kpc, reassuringly 



close to the tangent-point distance. Therefore if 
we adopt a mean value of (i=2.5kpc we would 
likely only need to attach a 12% uncertainty to it. 
Although their association in velocity and on the 
sky is strong circumstantial evidence, it is not cer- 
tain, however, that BYF73 is actually associated 
with the T] Car GMC complex or NGC 3324. In 
particular. Fig. [2] shows no evidence that the low- 
density molecular envelope of BYF73 has been at 
all disturbed by the vigorous star-formation activ- 
ity closer to r] Car or by the bubble of NGC 3324. 
Nevertheless, further evidence that the tangent- 
point distance is reasonable follows from analysis 
of the cm-continuum emission of the small HII re- 
gion adjacent to BYF73 (see Fi^ [Tlan d gMl). Us- 
ing the MGPS-2 ("Muruh v et al]l2007h and SGPS 
( Haver ko rn et al . 2006) fiux densities at 843 and 
1420 MHz of 62±5 and 85±11 mJy respectively, 
and assuming an electron temperature in the HII 
region Te = 7000 ± lOOOK (S haver e^il 19831) 
standard analysis iMezgei et al J Il967l : iBarnesI 
1985[ ) gives a distance-independent emission mea- 



sure EM= (1.4 ± 0.4) X lO^pccm"^, typical of 
compact HII regions fflabingfc lsraeil979h . Such 
HII regions have diameters ~0.1-lpc, bracketing 
that for BYF73 (from measurement of the Br-7 
nebula in Fig.[7l2, its FWHM = 0.25 pc), and so 
yielding a most likely location for it at the tangent 
point. 

In summary, various lines of reasoning make a 
good case for BYF73 lying close to the tangent- 
point distance for its longitude. We therefore as- 
sign a distance of 2.5±0 .3kpc, based on the di rect 
measurements listed bv lYonekura et all ( 20051 ). 
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3.3. Evidence for Gravitational Infall 

The dense molecular clump, centred at = 
(286°208,+0°169) and easily visible in the Mopra 
maps, has HCO+ sp ectral line p r ofiles that fit the 
canonical pattern of IZhou et al. I (|l993h indicating; 
gravitational infall onto a protostar (see Fig. 2]). 
For the optically thick HCO+ emission, this in- 
cludes a self-absorbed profile with predominantly 
stronger blue wings at most positions (the "blue 
asymmetry" or inverse P-Cygni profile, seen in 
panels a and c of Fig. [4]), together with more Gaus- 
sian line profiles for the optically thin transitions 
of H^^CO^, which are centred in velocity on the 
HCO+ self-absorption (panels h and d). Further, 
the J=4^3 lines (panels c and d) are brighter than 
the corresponding J=1^0 lines (panels a and 6), 
the self-absorption in the HCO+ is deeper in the 
J=1^0 than the J=4^3, the velocity difference 
Vhiue — Vred (sec bclow) between the blue and red 
peaks of the HCO+ lines is slightly greater in the 
J =1^0 than the J=4^3, and the blue and red 
peaks in the J=4^3 line are both slightly red- 
ward of the respective peaks in the J=1^0 line. 
All of these det ails are co mpletely cons istent with 
the lZhou et all (l993) and lMvers eraD (|l996) pic- 
ture of a dense core undergoing gravitational in- 
fall, where the velocity of the infall and the tem- 
perature both increase towards the centre, produc- 
ing the respective line profiles and ratios. 

However, the mass scale of the infall appears 
to be unusual. 



Mvers et al.l (|1996l ) developed a 



simple but useful two-layer model to evaluate ba- 
sic parameters from spectra of molecular cores 
which are undergoing_gravitational infall, while 
2 



De Vries fc MversI (|2QQ5l ) extended this analytic 
model and provide a general code for robustly de- 
termining these parameters and their uncertain- 
ties. Although these models were developed in 
the context of low- mass protostars, the results we 
derive here satisfy the assumptions made in their 
treatment of the radiative transfer. The key qual- 
ifications are that the infall speed not be much 
greater, nor much smaller, than the velocity dis- 
persion in the dense gas, which result we obtain 

I 1 1 1 

below. Here we use the M vers et al .' (1996; for- 
malism and the HCO+ line profiles to estimate the 
characteristic gas infall speed a nd motivate further 
discussion. In §3.41 we use the IPe Vries fc MversI 



From Myers et al. ( 1996l )'s eq. (9) and using the 
parameters as listed in Table [T] from the sample 
spectra in Figure [H we obtain Vin as shown also 
in Table [H and where the quoted errors are ob- 
tained by propagating the measured uncertainties 
through the formula. Continuing to follow My- 
ers et al., we need an estimate for the radius over 
which the infall profile is seen, in order to allow 
calculation of a kinematic mass infall rate. This 
profile is widespread in the HCO+ data, but its 
intensity drops only slowly into the background, 
showing no hard edge. To indicate a radius we 
consider the emission FWHMs (suitably decon- 
volved). In the HCO+ and R^^CO^J=1^0 data, 
the diameters Dfwhm = 120'' ± A" and 65'' ib 5" 
respectively, taking a geometric mean of the ma- 
jor and minor axes in each case. At 2.5 kpc these 
respectively give clump radii R = 0.73 ± 0.09 pc 
and 0.40 ± 0.05 pc, where we have now also added 
in quadrature the uncertainty due to the distance. 
Although the HCO+ infall profiles are clearly more 
widespread than the H^^CO^ radius, we conserva- 
tively take the latter as an optically thin tracer and 
therefore more representative of the true column 
distribution, understanding that this may in fact 
be a lower limit to the clump radius. This gives 



dt 



(2.9 ±1.5) X lO"^M0yr-^ (1) 



for BYF73' s mass infall rate, where jJimoi — 

2.30 

is the mean molecular mass in the gas and Tl/Qq^ IS 
the critical density for the J=1^0 transition (see 
^3.5p . This should be compared to the gravita- 
tional mass infall rate for the sel f-similar si ngular 
isothermal sphere (SIS) solution (|Shulll977l ) 

^ = ^^ 0.080 X 10-2 Moyr-\ (2) 

where instead of the sound spe ed c of Shu, we have 
substi tuted, as suggested by iBanerjee fc Pudritz 
(|2007l ). the supersonic velocity indicated by the 
H13C0+ linewidth a from Tabled] (see also 



(|2005l ) code to more rigorously evaluate the model 
fits. 



Even so, we see that for BYF73, Shu's solution 
cannot give us the observed infall rate. Instead, 
Banerjee fc Pudritz show that a magnetised core 
can collapse supersonically with an effective speed 
Mc, where M is the Mach number in the flow. 
For BYF73, then, the observed infall only requires 
collapse with M ~ 3. This Mach number and infall 
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Fig. 4. — Sample Mopra and ASTE spectra of BYF73. Panel a shows the HCO+ J=1^0 spectrum in 
black at the peak H^^CO^ position. The Hill5 model fit at this point (see §3.4p is shown in red, and the 
residual spectrum (data-model) is shown in green, offset IK below the HCO+ spectrum. The vertical blue 
line indicates the systemic velocity at Vlsr = — 19.7kms~^, also from the model fits. Panel b is the Mopra 
H-^^CO+ J=1^0 spectrum, also at the peak H^^CO+ position of BYF73, shown at the same and Vlsr 
scales as panel a. Panels c and d are ASTE J=4^3 spectra from the BYF73 peak position in HCO+ and 
H^^CO^, respectively, shown at the same velocity scale as panels a and 6, but on the scale. 



speed above are con sistent with (for example) the 
MH D simulations oflBaneriee fc Pudrit a (|2QQ7l ) or 
the iMcKee TanI ([2003') massive turbulent core 
model; however our mass infall rate is still more 
than an order of magnitude higher than in such 
models, mainly because of the large extent of the 
infall asymmetry in our maps. 

Mvers et"an ( 19961 ) suggested that, for their 



low-mass protostars, the agreement of the inferred 
and theoretical rates indicates the derived inward 
motions are consistent with gravitational infall. 
Under this interpretation BYF73 also gives a much 
larger infall rate than is typical of low-mass pro- 
tostars (~10~^ M0 yr~^, increasing to 10~^ to 
lO~^M0yr~^ during FU Orionis-type outbursts; 



Lada 19991 ). again stemming mainly from the 
parsec-scale extent of the asymmetric HCO+ pro- 
file, and also from the unusually large value for 
Vin- This mass infall rate is also larger than any 
seen so fa r in any similar massive star-forming; re - 
gion (e.g. iFuller et al.ll2005l : iBeltran et aTll2006h . 
Given the linear size of this region and the near- 
IR appearance of peculiar emission-line nebulos- 
ity at the centre of the clump, we suspect that 
the entire BYF73 cloud is undergoing a global 
gravitational collapse. Verification of this sugges- 
tion awaits additional supporting evidence includ- 
ing interferometric observations and more detailed 
modelling. However all of the Mopra spectral line 
maps of BYF73 (e.g. HCN, N2H+, etc.), as weh as 
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Table 1: 
the 



Sample infall fitting parameters for 



Myers etaJ] ( 19961 ) model at the peak 



-25 -20 -15 

Velocity (km/s) 



Fig. 5. — Sample Hill5 model spectrum (the same 
as Fig.[4]2 in red) with me asurements as shown fol- 
lowing Fig. 2c of lMvers"et al. (1996). 



the C S J=2^1 data reported bv lBronfman et aJ] 
(|l996l ). show similar emission distributions and/or 
line profiles, with differences as expected from the 
species' different relative abundance. This is not 
surprising considering that t hey all requ ire high 



densities {ricr ~ 10 cm~'^: iSpitzerl 1197 8) to be 
collisionally excited and thermalised to the gas ki- 
netic temperature, and so should reflect the same 
dynamical state as seen in the HCO+. 

3.4. Radiative Transfer Modelling 

iDe Vries fc Mvers (2005) compared a number 
of analytic radiative transfer models of infall in a 
low-mass dense core to a full Monte Carlo model. 
They found that their "Hill5" model gave the most 
accurate simulation of the Monte Carlo solution, 
and of all the analytic models they examined, was 
the most robust against various numerical and 
instrumental uncertainties. We have used their 
HILL5 code to analyse our Mopra HCO+ J=1^0 
data cube pixel-by-pixel, and present the results 
here. 

The parameters fitted by the Hill model to 
an infall spectrum include the peak line exci- 
tation temperature and optical depth, the in- 
trinsic (i.e., equivalent optically-thin) Vlsr and 
velocity dispersion, and the infall speed Vin- 
These fits, in ord e r to b e considered reliable by 
IPe Vries fc MversI ( 2005 ). should be to spectra 
with S/N > 30. Our Mopra HCO+J=1^0 cube 



Hi3CO+J=1^0 position 



Parameter 


J=1^0 


J=4^3 


Units 


Vblue 


-21.0 ±0.1^ 


-19.7 ±0.1^ 


kms""*^ 


Vred 


-18.0 ±o.r 


-17.8 ±0.r 


kms""*^ 


Td 


0.44 ± 0.02" 


3.6±0.r 


K 


Tbd 


1.10 ±0.04" 


2.6 ±0.2^ 


K 


Trd 


0.24 ± 0.04" 


0.9 ±0.2^ 


K 


<j 


1.5 ±0.3^ 


1.06 ±0.15^ 


kms""*^ 




0.86 ± 0.36 


0.34 ±0.12 


kms""*^ 



Footnotes: Measurements from (a) Fig.|4jz; (b) Gaussian 
fitting (not shown) to Fig.|4j?; (c) Fig.lJ]:; (d) Gaussian 
fitting (not shown) to Fig.|4ji. 



does not formally satisfy this requirement per 
0.11 km s~^ channel (peaking at S/N ~12), but 
since the fits are over a large velocity range (up 
to ~8kms~^ to zero-power) we suggest that the 
figure of merit should rather be the peak S/N 
in the integrated intensity map (Fig. [2]), which is 
^50. Put another way, by fitting 5 parameters 
to spectra with ~70 independent resolution ele- 
ments, the problem is more strongly constrained 
than the per-channel S/N would suggest. Our 
claim of reliability is bolstered by the solutions 
themselves, which show little statistical noise in 
the output parameters above an integrated inten- 
sity ~ 1.5Kkms~^, unless the optical depth is too 
low to give a good infall solution (i.e., everywhere 
except toward the SW edge of the clump). In such 
areas, the fitted T^^x and r especially are poorly 
constrained and not physical. 

In Figure [6] we show the parameter maps of the 
Hill5 model fits to the HCO+J=1->0 line. The 
lowest reliable excitation temperature (panel a) is 
^ 3-4 K around the northern and eastern perime- 
ter of the clump, and reaches a maximum ~6K 
near the peak emission. These values are clearly 
less than the T^.^ discussed below, but this is partly 
attributable to the data cube being on the T\ 
scale, which is 0.64 xT^^. The only effect this 
choice of temperature scale has on the model fits is 
on the scaling of the fitted T^.^ • Over the same ar- 
eas, the optical depth (panel 6) ranges from ~ 1.5 
to ~6. The areas where the infall solution is un- 
physical are clearly visible in these two panels as 
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(a) 



HillS model: T„ 




Galactic Longitude 



Fig. 6. — Pixel- by-pixel parameter fits of the HillS model to the Mopra HCO^ J=1^0 cube. All panels are 
overlaid with the same HCO+ integrated intensity contours from Fig.O (a) Peak excitation temperature. 
{b) Peak optical depth, (c) Intrinsic Vlsr- {d) Velocity difference between Fig.lSji and panel c in this figure, 
(e) Intrinsic velocity dispersion. (/) Infall speed. 



the red or black pixels, predominantly to the south 
and/or west of the clump. In the brighter ar- 
eas where we believe the solutions are reasonable, 
we note that the highest optical depth lies NE 
of the brightest HCO+ emission, precisely where 
the brightest H^^CO^ ridge lies. Furthermore, the 
HCO+/H^^CO+ line ratios (ranging from ^6-10; 
see Fig.[2j) are entirely consistent with these opti- 
cal depths and normal abundance ratios of these 
molecules. This is actually remarkable since the 
Hill5 model only fits parameters to the HCO+ 
cube: the fact that the H^^CO^ data are consis- 
tent with the model results gives us further con- 
fidence in the HillS solutions. Likewise, the high- 
est Tex is to the SW of the HCO+ peak, approx- 
imately facing the HII region (cf. § §3.2|3.6p and 
entirely consistent with that geometry. 

The various velocity parameters of the HillS 
models are also remarkably well-behaved. The 
systemic Vlsr map (panel c) looks grossly sim- 
ilar to the moment- 1 map (Fig.[3^), but in fact is 
slightly redshifted where the infall profile is most 
prominent. This is made clear in panel d which 



shows the velocity difference (moment- l)-(VL5i?). 
This quantity should be close to zero in most 
places, but skewed to negative values where the 
infall is strong and the moment- 1 values reflect the 
blue asymmetry of the spectra. Indeed the colours 
in panel d show exactly this: away from the robust 
infall solutions, the average colour is orange corre- 
sponding to a mean difference ~Okms~^. Where 
panels a and b have good solutions, the mean ve- 
locity difference is consistently ~ — 0.5kms~^, in- 
dicating the extent of the spectral asymmetry. 

The last two panels of Figure [6] show the 
velocity dispersion (panel e) and infall speed 
(panel /). In the area of good fits, the for- 
mer is ~ 1.0±0.2kms~^, while the latter is 
~ 1.0±0.4kms~^. Once again, we see these 
HCO^-derived dispersions are consistent with the 
actual linewidth measurements of the H^^CO^ 

i 33.1D . Furthermore, as predicted bv'Pe Vries fc MyersI 
2005) the Hills solutions do indeed scale to 
higher-mass regions than they examined, since for 
BYF73 the criterion that the intrinsic dispersion 
is comparable to the infall speed is satisfied. 
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In panel / we note an interesting structure in 
the velocity field of the infall. Along the long axis 
of the clump (i.e., to the NW and SE), the in- 
fall speeds are consistently lower than lkms~^, 
whereas across the short axis (to the NE and 
SW) the infall speeds are consistently higher than 
lkms~^. It is tempting to interpret this pattern 
as due to a partially rotationally-supported oblate 
clump. In this scenario, the infall is somewhat cen- 
trifugally hindered in the equatorial plane (which 
is roughly parallel to the major axis of the emis- 
sion) by a rotational speed which may be ~0.5- 
lkms~^, but is unimpeded along the presumed 
rotational axis (roughly parallel to the minor axis). 

Equally, we note that the Vlsr field in panel 
c is purely kinematic, since any radiative transfer 
effects would have been filtered by the model into 
just Tea; and r (panels a and b). Thus panel c 
may be a better indicator of rotation in BYF73, 
with the rotation axis being roughly aligned with 
the emission's major axis instead, and suggesting 
a more prolate geometry for BYF73. However we 
also note that the most redshifted portion of panel 
c, centred near (286° 195, 0.°L57), has already been 
attributed to a separate, non-infalling component 
in the data cube ( §3.ip . Both of these rotation in- 
terpretations are thus quite speculative: the noise 
in both the data and the model may dominate the 
features we are trying to interpret, and these al- 
ternatives really need to be explored with higher 
resolution and greater sensitivity data in order to 
discern between them. 

An important feature of the Hill5 Vin treat- 
ment is that it is consistent with the val ue ob - 

( 1996h 



based on the iMvers et al 



tained in 

work, w hich was for a two- l ayer r adiative transfer 
model. IPe Vries fc MversI (|2QQ5[ ) similarly found 
that their two-layer models often gave solutions 
for Vin which were mostly consistent with the Hill 
models, however the Hill5 model was the most ro- 
bust to errors. In the calculations below we use 
the Hills result Vin = 1.0±0.4 kms-^ For the 
mass infall rate from §3.3[ we now have a some- 
what larger value 



dMk 
dt 



(3.4 ± 1.7) X lO-2M0yr-\ (3) 



recalling that by evaluating this with the H^^CO^ 
radius we are likely obtaining a lower limit to the 
global infall rate. We conclude that the radiative 



transfer modelling of the HCO+ data gives results 
which are surprising but highly self-consistent, and 
consistent with other features of our data. 

3.5. Clump Mass 

Despite the satisfactory results of the mod- 
elling, in order to make a strong case for the 
formation of a massive cluster, we also need to 
establish that the molecular clump has sufficient 
mass to qualify for this status, and that other pos- 
sibilities for interpreting our data are discounted. 
Since the gas density will probably be at least 
the HC0^J=1-^Q transition's critical density 
(jHaese fc Woodsl Il979l : iBarnes fc Crutcherl [T99Q[ ) 
where the bright molecular emission is seen, the 
cloud mass is given approximately by 



M > iiimoimHncr{7r/ln2f^^R^ 



(4) 



1.0 X 10^ 



6.4 X 10^ Mr. 



3 X lO^cm-3 J V0.40pc 



R 



3 X lO^cm-3 J V0.73pc 



R 



using the volume for a 3D Gaussian. Here we give 
two values for the mass based on which size we 
take for the HCO+-emitting region. With ASTE's 
detection of the J=4^3 line, even higher-density 
gas (~10^ cm~^) must exist in the clump, and if 
widespread would give a much higher mass esti- 
mate. Therefore the first value for the mass is 
almost certainly a lower limit. However eq. (4) as- 
sumes that the dense gas giving rise to the emis- 
sion fills our beam, whereas the filling factor / is 
unknown and possibly <Cl; this may indeed be the 
case in the outer envelope of the clump, thus the 
second value is probably an upper limit. 

A formally more rigorous, but not necessarily 
more precise, mass estimate is made (and we ob- 
tain an estimate for / as well) if we calculate the 
HCO+ column density first. We use the full ex- 
pression without assumptions about optical depth 
or approximations to the stimulated emission cor- 
rectio n in the denominator (e.g.. jRohlfs fc Wilson 



2006D . Assuming LTE applies and with quantities 
in cgs units, we obtain a column density for each 
line of sight from 



7V(HC0+) 



3h Q(Te,)e^-/^^- 



rdV , 
(5a) 
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where Q is the partition function for HCO+ at 
the excitation temperature Tga^, Eu is the energy 
level of the upper state Ju of the transition, /i^ is 
the molecule's electric dipole moment, and the line 
optical depth r (peak value ~6 from the previous 
section) is integrated over the velocity, here taken 
over the range -23.2 to -16.6 kms~^ (as in Fig. [2]). 

Determining the excitation temperature is a 
little more complicated, however. Faundez et 
al (2004) derive = 30 K for the continuum 
dust emission from BYF73, but found it necessary 
to fit two temperature components to the spec- 
tral energy distributions (SEDs) of most of their 
sources. They do not give an explicit value for 
the warm component in BYF73, but their aver- 
age warm component has Td ^ 140 K. From the 
HCO+J=4^3/1^0 brightness ratio (2.20±0.14 
from Fig.m suitably corrected for the beam ef- 
ficiencies) at the peak of BYF73, we fit a 
= 125±26K for the dense gas. But without a 
spatially-resolved HCO+ J=4^3 map, we are lim- 
ited to saying that the gas T^x probably takes a 
range of values from 30-125 K. Because HCO+ is 
a linear molecule, its partition function is stra ight- 
forward to calculate ( Rohlfs fc Wilson 

^M). At 

these temperatures Q ~ 14 — 59, giving 



total cloud mass 

Mlte = -TpdJ-moirnHj-r-^r 
A InZ 

- (1.2- 17) X 10'' Me 



(7) 



{R/OAOpcf 
(X/10-9) ■ 



iV(HCO+) 



4.84 X 10 



1 — ^-hu/kTex 
15 ^^-2 



/ 



rdV 



(0.92-13) X 10'^ cm 



where the velocity is in kms~^, and we have taken 
a Gaussian line profile, with dispersion a as be- 
fore, for the integral. Combining the column den- 
sity with the size measurement (assuming that the 
physical depth of the source is comparable to the 
projected size) gives a central density estimate 



ln2 N 



The lower limits for both eqs. (6) & (7) are prob- 
ably too low, since they don't include the contri- 
bution to the mass from the warmer component; 
moreover we have used the smaller H^^CO^ ra- 
dius in both. With the larger radius, eq. (6) gives 
a density ~1.9x 10^ cm~^, lower than before but 
still near the critical density, and eq. (7) a mass 
~4.1xlO^M0. Both low-temperature mass esti- 
mates are close to the values in eq. (4). This is 
surprising given the approximate nature of these 
calculations, but reassuring. We conclude that an 
intermediate value, M ~ 2.OxlO^M0, is probably 
reasonable for BYF73. 

However the upper limits from eqs. (5-7), 
based on the higher excitation temperature be- 
ing widespread, are certainly too high, since 
it is unlikely that such a warm temperature 
would be typical of the whole parsec-wide dense 
clump. One would need to compare a map of the 
HCO+J=4^3 emission at a resolution at least 
as good as our HCO+J=1^0 map, in order to 
cm~^ derive a reliable map of T^x across the source and 
address how much larger the clump's mass might 
(5b) be due to this warmer gas. 

In eq. (6-7) we have used an abundance X = 
10~^ for HCO+ relative to H2, which is a strong 
upper limit from so me recent models of massive 
core chemistry (e.g.. iGarrod et al]l2008l ). These 
models show X^qq+ is a strong function of time, 
and is not necessarily the main charge carrier in 
such regions. Thus Xhco+ massive cores may 
be an order of magnitude lower than a more typical 
value ^ 10~^ in l ow-mas s cores (e .g., Lor en et al. 



(3.5 



RX 

51) X lO^cm-3 



(6) 



IQQdlCaselh et a l. 2002: Lee et al. 2003) an d used 



(i?/0.40pc)(X/10-9) 



over the same temperature range, which shows 
that the central density in BYF73 almost certainly 
exceeds the critical density for thermalising the 
HCO+J=1-^0 line, and vindicates this assump- 
tion in eq. (4). Similarly integrating the column 
density (eq. 5b) over the emission region yields a 



here. On the other hand, Zinchenko et all (|2009l ) 
obtain Xhco+ ^ 2.3-12 x 10~^ cm~^ from obser- 
vations of a sample of massive clumps, although 
they cautioned that their values are probably over- 
estimates. This means that parameters derived 
here that depend on X are quite uncertain. Nev- 
ertheless, eqs. (4) and (7) suggest that / may be 
close to unity, and that M is large. 

Indeed BYF73 seems to be quite extreme in this 
regard too. For example, the mass surface density 
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corresponding to the column density from eq. (5b) 
is S = {N / X){iinioi^H) ~ 35kgm~^, which is 
near the largest value of t he m assive Gal actic clus- 
ters considered by McK ee fc TanI (|2QQ3h . There- 
fore BYF73 is interesting as a likely environment 
in which massive protostellar cores may form, and 
then form massive stars. 

Are there alternatives for the dynamical state 
of this clump besides gravitational infall? To an- 
swer this, we evaluate a number of terms from 
the Virial Theorem. If the linewidths seen in 
the HCO+ (~2kms~^ relative to the line centre, 
counting emission out to the half-power level) were 
due to rotational support against self-gravity (an 
interpretation we do not favour due to the self- 
consistency of the infall modelling, and the effec- 
tive limit of ~ 1 kms~^ to any rotational contribu- 
tion to the spectral lines), then 



Mrot = V^R/G 

- 370 Mo 



^a) 



2 kms 



R 



0.40 pc 



However thermal and magnetic pressure must also 
contribute to the support of the cloud; the corre- 
sponding virial relations give 



Mth = ^kTexR/rriH^G 

Tex 



(8b) 



240 M, 







R 



125 Ky V0-40pc 



and (in cgs units only) 



mag 



1/2 



(55^i?718G) 



(8c) 



940 Mo 







600 /iCy V0-40pc 



R 



where we have taken an appropriate value from 
studies in similar regions for the magnetic field, 
up to twice the typical Zeeman-derive d Bias at 
a de nsity 3xl0^cm~^ (see Fig. 1 of ICrutcherl 
19991). Such a large va lue is further supported by 
Falgarone et al.l |2008 l) who obtained 5 (median) 
= 560/iG from CN observations of clouds with 
a mean density 4xl0^cm~^. These terms, even 
in combination (~1550 M©), are much too small 
to provide the necessary support against gravity, 
unless (for example) the magnetic field strength 
were at least ten times the value assumed here, 
and/or we take linewidths out to the zero-power 



level (±4 kms~^). While such values for rotation 
and the magnetic field are not entirely ruled out 
as a means of supporting BYF73 against collapse, 
they would be quite extreme. We infer that virial 
equilibrium does not apply in this case, despite the 
HCO+ abundance and T^x uncertainties. 

There is also the possibility that the velocity 
pattern in BYF73 represents a massive outflow 
rather than infall. Besides the detailed spectro- 
scopic arguments for infall, we discount the out- 
flow interpretation since maps of the HCO+ line 
wings (not shown here) do not reveal any particu- 
lar geometric pattern, such as a bipolar separation 
of the line wings. Nevertheless, sensitive ^^CO ob- 
servations should be made of BYF73, since they 
might be better able to find any outflow, if present. 

The conclusion that BYF73 is indeed a mas- 
sive dense clump undergoing contraction at least 
(if not collapse) seems fairly reliable, the strongest 
evidence being the line profiles, the mass calcula- 
tions, and the IR appearance (see §3.6p . We can 
compare our mass estimate for BYF73 with oth- 
ers in the literat ure. From Nant e n CO mapping 
and IRAS fluxes, lYonekura et al.l (|2005l ) obtained 
LTE and virial cloud masses of 1900 and 3600 M© 
(resp.) and a luminosity 3.OxlO^L0. Our LTE 
mass estimates are significantly larger than theirs, 
likely due to our use of a tracer of denser gas, 
but our total virial mass is smaller than theirs, as 
might be expected with a smaller observed linear 
size and linewidth. Our LTE estimates would be 
smaller if we assumed a smaller effective density 
for HCO+ and/or a larger IICO+ abundance. Ei- 
ther of these might reduce our best estimate above 
by a factor of 3 or so , to <10^ M^;^, br ings:ing it more 
into line with the l Yonekura et al. I 12005.) num- 
ber, although we do not favour th is value. From 
1.2mm mapping and SED fitting iFaundez et al.l 
(|2004l ) obtained a clump mass of 470 M© , density 
= 1.4xl0^cm~^, dust temperature 30 K, and 
luminosity 1.9xlO^L0. Their mass value seems 
quite low, but could be as much as five times 
higher with a lower assumed dust opacity, as they 
point out. Thi s would bring it into c loser agree- 



ment with the Yonekura et al. ( 20051 ) mass esti- 



mate, and suggests that such lower dust opacities 
may be required to explain the higher molecular 
masses. At any rate, our IICO+ data strongly sug- 
gest the presence of a large amount of dense gas 
that may not be fully sensed in CO lines or in the 
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Fig. 7. — (a) RGB-pseudocolour image of BYF73 in i^-band spectral lines. Here Br-7 is shown as red, and 
H2 is shown as green {v=1^0) & blue {v=2^1). Contours are overlaid from Mopra HCO^ (gold) and 
H^^CO+ (red) integrated intensities (levels as in Figs.[3]&[2l resp.). (b) Same contour levels as a, but now 
coloured blue and red^ and overlaid on the i^-continuum image. The stars labelled in magenta are the same 
as those similarly labelled in Fig.lHK- Recall from Fig. 1 that at a distance of 2.5 kpc, the scale is 40" = 
0.485 pc or 0°.02 = 0.873 pc or Ipc = 0°.0229 = 82^:5. 



mm-continuum. 



3.6. Infrared Features 



Our i^-band imaging (Fig.[7j) has much higher 
angular resolution than our mm data, or archival 
centimeter- wave (cm) and far-infrared (FIR) data, 
and shows some striking structures and correla- 
tions. Near the molecular clump, there is a com- 
pact HII region (visible as a Br-7 emission nebula) 
and IR cluster (visible also in the i^-continuum 
image). The Br-7 is exactly coincident with a 
centimetre-continuum point source from both the 
Molonglo Galactic Plane Survey-2 (the MGPS-2 
has a sim ilar beamsize, 45" x 53'', to our Mo- 
pra data; iMurphv et al. 1 l2007h and the slightly 
l ower - resolution South ern Galactic Plane Survey 
(|Haverkorn et al.ll2006l ). Such features have been 



seen before around similar massive star-formin i 
dense clumps, e.ff^GC 202 £jBarnes et al.lllQs" 
or AFGL 5179 (j^ et allbooi ). but the example 
of BYF73 is interesting in the rather clean sepa- 
ration of the ionised and molecular components, 
and the distinct "cocooning" of the excited H2 
emission around the very symmetric Br-7 cm- 
continuum. Indeed, this is actually reminiscent 
of planetary nebulae (e.g.. Fig. 9 of iRvder et al 



19981 ) or the classic picture of a Stromgren sphere. 
From the pseudocolour composite image in Fig.[7K, 
we see that the shell of excited H2 appears to 
be traced much better by the v=2^1 than the 
'^=1^0 emission. This is surprising since a [1- 
0]/[2-l] ratio less than unity would be at odds 
with our understanding of H2 excitation. Instead, 
this ratio is likely either an artifact of differential 
reddening between the various filters used, or due 
to non-photometric imaging conditions, or both. 

To confirm this we obtained a iC-band long- 
slit spectrum aligned between the mm and Br-7 
peaks (Fig. [8]), which shows that the [1-0]/ [2-1] ra- 
tio is actually 2.14±0.10 at the molecular- ionised 
interface, typical of photodissocia tion regions 



(PDRs) or shock-excited jets (e.g. , lAllers et al 



2OO5I : ICaratti o Garatti et al.l l2006h . From the 



measured S'(l) [l-0]/[2-l] ratio and the tabulation 
of T. Geballe (1995jj, the gas kinetic temperature 
at the PDR interface is constrained to be >4000K; 
including the S{0) and S'(3) lines also visible in 
Figure [8] (with respective ratios 5'(1) 1-0 to these 
lines of 1.65±0.07 and 4.5±0.3) suggests a tem- 



s Quoted by T.Kerr (2004), www.jach.hawaii.edu/UKIRT/ 
astronomy /cahb /spec_cal /h2_s. ht ml 
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Fig. 8. — (left) Near-IR long-slit spectrum across BYF73, where the slit was oriented along a line joining 
the peak of the HII/Br-7 emission (to the NW of Fig. [71 at the top of this figure) and the peak of the HCO+ 
emission (to the SE of Fig. [71 at the bottom of this figure). The spectral lines are as indicated, and the labels 
A, B, and C refer to stars similarly labelled in Fig. [3 (right) Sample spectra from panel (a) on a relative 
fiux scale to indicate line ratios. The upper spectrum (labelled "SE") corresponds approximately to pixel 
row 340 in panel (a), while the lower spectrum (labelled "NW") is near pixel row 200. 



perature as high as ~5000K. This is comparable 
to, but somewhat higher than, H2 t emperatures 



seen in other star format ion PDRs (lAllers et al 



2005 ) or low- mass H2 jets (|Caratti o Garatti et al 
2006h . approaching the typical Tg - 7000 K for 



Galact ic HII regions at th is galactocentric radius 
(8 kpc: [shaver et al]|l983[ ), and is perhaps indica- 
tive of the relative youth of the HII region in 
BYF73, and/or the strength of the shock excita- 
tion from the young stars in the HII region. 

Another remarkable feature of the IR imaging 
is the apparent deficit in line emission from the ex- 
act peak of the mm-molecular emission (within the 
pointing uncertainty). Unfortunately a spectrum 
for this position is not available, since it coincides 
with emission in the reference beam (the horizon- 
tal black line in Fig. [HI at pixel coordinate 420). So 
while we cannot obtain any IR line ratios here with 
the current data, it appears as if the H2 v=2^ 1 
and Br-7 are both seen in absorption at this peak, 
creating an apparent "absorption nebula". This 
nebula can be seen in Figure [T^i as a green patch 
to the left (Galactic east) of the HII region, since 
there the blue (H2 v=2^1) and red (Br-7) ap- 
pear more strongly "absorbed", while the green 
(H2 v=1^0) is only weakly "absorbed". While 



much of this appearance may be due to the con- 
struction of the RGB image, at the very least it is 
likely that there is either unusual, highly localised 
IR emission/absorption at the molecular peak, or 
that the deeply embedded stars at that position 
have very unusual IR colours. Moreover, this po- 
sitional coincidence is highly suggestive. Further 
east of the peak of the absorption nebula, there 
seems to be a weaker, comma-shaped extension, 
as well as a highly reddened cluster of stars; this 
shape is also seen in some of the HCO^ channel 
maps. Furthermore, this extension is also aligned 
with the supposed equatorial plane of the oblate 
clump scenario from §3.41 

There is evidence in other archival data for 
an unusual source at the molecular peak. A 
3-colour (i.e., JHK) 2MASS image shows that 
the HII region exhibits some moderate redden- 
ing, but that the knot of K-band emission visi- 
ble in Figure [7f) at the molecular peak is virtu- 
ally invisible at shorter wavelengths, confirming 
its highly embedded nature. A similar 3-colour 
MSX image (i.e., 8,15,21 /im) further shows that 
this embedded source dominates the luminosity at 
MIR and longer wavelengths. We therefore have 
the rather unusual situation that, while a more 
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evolved source is adjacent to our molecular clump, 
the apparently less evolved source (s) within the 
clump are more luminous than the revealed ex- 
citing stars. Indeed, the luminosity of the deeply 
embedded IR source (s) is at least partly derived 
from the release of gravitational potential energy. 
Taking the mass inflow rate from §3.4[ which is 
itself probably a lower limit as described there. 



■^grav 



GMM 
R 

1200 



0- 



(9) 

(M/20,000MQ)(M/0.034MQyr-^) 
(i?/0.40pc) ' 



This result, that Lgrav ~ ^%LhoU i^^iy be even 
larger if higher-resolution observations of the in- 
fall and central MIR/NIR sources reveal the inflow 
continues deeper into the central regions. 

If the NIR "absorption" at this position were 
real and not an imaging artifact, then it would im- 
ply very large columns of gas, ~10^^cm~^. The 
HCO^ self-absorption and spatial distribution re- 
quire the same condition, and so the implication 
of very high column density at the molecular peak 
would seem to be strong. This is further supported 
when we note that the continuum emission of the 
three stars closest to the "absorption nebula" po- 
sition in Figure [8^1 (i.e., the bright horizontal lines 
at pixels 400, 435, and 460) show a very strong at- 
tenuation at the shorter i^T-band wavelengths, pre- 
sumably due to severe reddening in the molecular 
clump. Such reddening is not apparently affecting 
the stars in the HII region (e.g., those labelled A, 
B, & C in Fig. [8]) to the same degree. 

The i^-continuum image reveals details of clus- 
tering in BYF73. Compared to the surrounding 
sky away from any HCO+ emission, within the HII 
region there is clearly an overabundance of stars. 
In addition, the bright H2 nebulosity immediately 
to the east of the HII region contains an even more 
compact clustering of brighter stars, and there is 
another tight grouping around the molecular peak. 
To the north and south of the molecular peak, the 
star density is actually lower than the surround- 
ings, suggesting that here the dust column den- 
sity is still so high that background stars are be- 
ing extinguished at 2/im. It is clear that many 
young stars, some massive enough to form an HII 
region, have already formed in BYF73, and that 
further star formation appears to be proceeding 
vigorously to the east of this HII region. 



3.7. Theoretical Considerations 

We note that the typical projected nearest- 
neighbour separations of stars in these IR groups 
(i.e., at the molecular and Br-7 peaks, and molecular- 
ionised interface), ^2" or 5000 AU by inspection 
of Figure [7)3, is less than the Jeans length 



Rjea 



kin 



GifJ-moirriHYnH^ 
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3 X 10^ cm- 



(10) 



-0.5 



in the dense clump, and although higher densities, 
especially near the centre, may make these scales 
more commensurate, our estimate for eq. (10) is 
probably a strong lower limit considering that 
there is warmer gas in the clump ( ^3.5p . and 
that the Imm-der ived density (1.4x 10^ cm~^, 
Faundez et al ] l2004h is lower than that used above. 



This disparity is typical of massive clusters (e.g. 
Churchwellll2QQ2l ) and is a well-known feature of 
such regions that models must reproduce. Cur- 
rently, theories attempt to model this structure 
using either co mpetitive processes (such as co- 
alescence, e.g. Bonnell et al. l2003[ ) or a scaled 
up accretion disk / turbulent core scenario (e.g., 
McKee fc Tanll200 3h. BYF73 promises to be a use- 



ful test case in this debate, but as suggested by the 
IR imagery, will require mm-interferometric obser- 
vations that approach the IR resolution. At this 
level (~1'' or better) we begin to match (at the dis- 
tance of BYF73) the spatial resolution (^ .01 pc ) 
in the simulations of iBaneriee fc Pudrit j (|2007h . 
High-resolution maps of the gas velocity field and 
linewidth will then help to discriminate between 
the competing theories. 

From the current mm data we can say that 
the line emission pattern and derived velocity field 
in BYF73 are consistent with the detailed MHD 
simulations of Banerjee fc Pudrit j (|2007) a nd ra - 
diative transfer treatment of Zhou et al.l (|l993h 
for protostars, as well as with the treatment of 
^cKee fc Tan (2003). However in this case (a) 
the mass infall rate and mass fc size scales are 
much larger than in any of these models, (b) there 
are multiple protostellar objects within the col- 
lapse zone, rather than a single, more massive 
one, and (c) the canonical spectral energy distr ibu- 
tion of low-mass Class protostars (Andre et al.l 
|2000[ ) has the flux dropping to undetectable lev- 
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els at wavelengths shortward of 10/im, although 
this is under the assumption of spherical symme- 
try. In BYF73 there are a number of near-IR 
sources visible at the centre of the infalling clump, 
so it is likely that orientation effects play a role in 
the emergent SED for massive protostars, and/or 
that the SED evolution is different in the massive 
protostar case. Again, higher-resolution mm- and 
FIR-continuum images of the cluster sources will 
help delineate SED evolution in these massive pro- 
tostars. 

Furthermore, this ongoing star formation is 
happening within a large-scale infall region (~1 pc), 
and to our knowledge this is the first time that 
such a coincidence of phenomena has been seen. 
With so much gas still infalling, it is possible 
that BYF73 could form many more stars before 
the supply of material is consumed. Dividing 
the clump's mass by the infall rate gives a max- 
imum lifetime t ^ 6 x 10^ yr for the supply of 
raw material for new stars, although if a proto- 
star massive enough to develop its own HII re- 
gion ionises the gas and arrests the infall, the 
cluster's formation may be complete in much less 
time. This is very long compared to a dynamical 
timescale, tff = 7r(i?^/8GM) ^2 ^ 30,000yr, and 
ten ds to support lon ger-timescale models such as 
the lTan et all (|2QQ6h "Equilibrium Cluster Forma- 
tion" model, rather than (e.g .) the r a pid star clus - 
ter formation models of Elmegre^ (120001 . 2007 ). 
However the most embedded stars in the IR "ab- 
sorption nebula" seem to be arranged in a fil- 
amentary geometry; being non-symmetric, such 
distributions tend to support a rapid formation 
scenario, however this appearance may be affected 
by the strong extinction in the area. 

In summary, we claim that the molecular 
and IR observations of BYF73 indicate the ex- 
istence of a dense clump undergoing global gravi- 
tational infall, simila r to the case for NGC 2264C 
(jPer etto et al.l 2006 ). We further suggest that, 
on the basis of the size, mass, luminosity, rate of 
infall, and near-IR appearance, BYF73 is in the 
process of forming a massive protocluster. The 
global mass infall rate as determined from the 
Mopra mm observations is very high even for a 
massive protostar, (3.4 ± 1.7) x lO~^M0yr~^ or 
more. To our knowledge, the upper end of this 
range would be unprecedented, if confirmed. 



4. Conclusions 

From Mopra and ASTE HCO+ observations, 
the Galactic source G286.21+0.17 (which we also 
call BYF73 from the CHaMP survey master list) 
has been found to be a massive dense molecular 
clump exhibiting clear signs of gravitational in- 
fall. The size and scale of this infall, dM/dt ~ 
3.4 X lO~^M0yr~^ over ~lpc, is either a record 
or close to it, and may indicate the global for- 
mation of a massive protocluster. A AT near-IR 
imaging confirms the existence of unusual spec- 
tral signatures and a deeply embedded cluster of 
stars in the infall zone, as well as an adjacent com- 
pact HII region and young star cluster. Higher- 
resolution mm- wave and FIR/MIR observations of 
this source are encouraged, since it appears to be 
an exemplary test case for confronting competing 
theories of massive star formation. 
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